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We study the physical properties of a half-metallic ferromagnet| superconductor (HM|S) bilayer, allowing for 
an arbitrary bulk pairing symmetry of the superconductor and spin-dependent processes at the interface. In 
particular, we study how the possibility of unconventional pairing such as p- and d-wave and a spin-active 
interface influence the (i) conductance spectra, (ii) proximity effect, and (Hi) local density of states of such a 
bilayer. Our calculation is done both analytically and numerically in the ballistic limit, using both a continuum- 
and lattice-model. It is found that the spin-dependent phase-shifts occuring at the HM|S interface seriously 
influence all of the aforementioned phenomena. We explain our results in terms of Andreev reflection in the 
presence of a spin-active interface, allowing for both spin-filtering and spin-mixing processes. We demonstrate 
how the surface-bound states induced by the anisotropy of the superconducting order parameter at the HM|S 
interface are highly sensitive to these spin-dependent processes. Our results can be directly tested experimentally 
using STM-measurements and/or point-contact spectroscopy. 



PACS numbers: 

I. INTRODUCTION 

In recent years, the physics of composite 
superconductor! ferromagnet systems has been subject to 
intense investigations. Apart from a wealth of interesting 
effects to explore from a fundamental physics point of view, 
it is also hoped that the interplay between the dissipation- 
less current flow in superconductors combined wifli tiie 
spin-polarization in ferromagnets will lead the way to new 
applications in low-temperature nanotechnology. 

The mutual influence of superconducting and ferromag- 
netic elements in heterostructures has a long history, see 
and references therein. While the basic constituent in 
a superconducting condensate is a spin-singlet Cooper pair 
in the usual Bardeen-Cooper-Schrieffer^ paradigm, the su- 
perconducting coiTelations are strongly altered when placed 
in close proximity to a ferromagnetic system, which spon- 
taneously breaks time-reversal symmetry, whenever transla- 
tional symmetry or time-reversal symmetry is broken. Cooper 
pairs with unconventional pairing correlations are formed in 
general'* ^. Such pairing correlations are unconventional in 
the sense that they differ from the conventional spin-singlet 
Cooper pairs, and they may exhibit for instance a spin-triplet 
symmetry or an odd-frequency symmetry. The study of the 
proximity effect in superconductorjferromagnet heterostruc- 
tures has received a lot of attention in recent years (see, e.g., 
Refs.EHlTll. 

In the extreme ferromagnetic limit of a half-metal, where 
the spin-polarization is close to 100%, one would naively ex- 
pect proximity-induced superconducting correlations to be de- 
stroyed due to the large exchange field in the ferromagnet. 
However, quite surprisingly at the time, Keizer et at. founcP 
that a supercurrent could flow between two conventional s- 
wave superconductors separated by a half-metallic layer of 
considerable size (^ /jm). This finding prompted several au- 
thors to investigate the underlying physics that permitted the 
superconducting correlations to survive over a long distance 



in a fully polarized ferromag nePSl. 

Prior to the experimental finding in Ref.''^, the concepts of 
spin-mixing and spin-flip processes were drawn upon in Ref.^ 
in order to explain how a supercurrent could be generated 
and sustained in an i-wave/half-metalAs-wave junction. The 
scattering of quasiparticles at the interface may in general be 
spin-dependent in the presence of magnetic parts of the sys- 
tem, which is the case for a superconductor/feiTomagnet junc- 
tion. This renders the transmission probabilities for spin-^ 
and spin-j^ particles different not only in magnitude, but also 
through the phases they pick up upon scattering at the inter- 
face. This gives rise to a so-called spin-mixing at the inter- 
face, which allows the singlet amplitude to be converted into a 
5- = triplet component, since scattered electrons with oppo- 
site spins experience different phase shifts at the interface. As 
a result, the superconducting correlations become a superposi- 
tion of both singlet and ~ triplet pairing. It is convenient 
for later use to briefly recapitulate here how this happenJ^] 
Consider a singlet correlation function in the superconductor: 

w) = n)k\i)-k-\i)kn)-k- (1) 

Upon scattering at the interface, the spins acquire different 
phase shifts 

lt)-i = e'*'lt)^, \i)-k=^''m)k- (2) 
This transforms Eq. ([T} into 

lv) = -cos(Ae)(|t)^i;>-i-i;>^lt>-A-) 

-isin(Ae)(|t>A.|;>-^ + i;>ilt>-^). (3) 

Here, A0 = 0| — 9|. The spin-dependent phase-shifts at the 
interface induce a triplet component which contributes to the 
total wavefunction \\\t) as long as AG ^ 0. 

However, it is also necessary to generate an equal-spin pair- 
ing S'j = ±1 components in order to sustain the long-range 
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triplet correlations. This demands spin-flip scattering pro- 
cesses of the type — > \l.)k and close to the inter- 
face. Such processes are unavoidably present for instance in 
the case where there are local inhomogeneities of the magnetic 
moment near the interface. The combination of spin-mixing 
and spin-flip processes then explain how the spin-singlet s- 
wave component of the bulk superconductor may be converted 
into a long-range spin-triplet component that is able to survive 
the large exchange field in the half-metallic region. 

The above discussion underlines the crucial importance of 
treating the interface properties correctly, and specifically tak- 
ing into account the spin-dependent phase-shifts that may oc- 
cur for the particles participating in the scattering processes^Sl 
In addition, the presence of mixed-parity pairing correlations 
in a S|HM structure should be linked to the spin-active nature 
of the interface. Previous literature has considered only the 
proximity effect between conventional s-wave superconduc- 
tors and half-metallic ferromagnets"^. In the present work, our 
aim is to investigate the interplay between the spin-dependent 
interface properties and unconventional pairing symmetries in 
the bulk superconductor with regard to the (i) conductance 
spectra, (ii) proximity effect, and (Hi) local density of states 
of such a bilayer. These quantities are directly accessible 
in experiments via STM-measurements and/or point-contact 
spectroscopy. In particular, by allowing for an unconventional 
pairing symmetry in the superconductor, such as p-wave or 



c/-wave, we may investigate the interplay between Andreev- 
bound surface states'*"^' and half-metallicity. 

We organize this work as follows. In Sec. |ll] we present 
the theoretical formulation used in this work, namely the 
Bogoliubov-de Gennes formalism 



III and IV 



we 



In Sec 

present and discuss our results for the conductance and prox- 
imity effect/DOS, respectively. Finally, we give our conclu- 
sions in Sec. |V] We will use boldface notation for 3-vectors 
for 4x4 matrices, and ^ for 2 x 2 matrices. 



II. THEORY 

In order to calculate the conductance of the S/HM junc- 
tion, we apply a modified Blonder- Tinkham-Klapwijk (BTK) 
theory which takes into account both an arbitrary pairing sym- 
metry of the superconductor as well as spin-mixing at the in- 
terface. Specifically, we consider the situation as shown in 
Fig. [T| where the region near the interface is allowed to have 
misaligned magnetic moments as compared to the bulk of 
the half-metallic ferromagnet. Our starting point is the BdG- 
equation 



H^> = e^' 



(4) 



in the half-metallic and superconducting region. We find that 



H = 



/Ho-h,@{-x)+V^d(x) 


A(e)*©(x) 



{V,-lVy)b{x) 

Ho + h,&{-x) + Vid{x) 

CA(e)*0W 







CA(0)0(^) 

-HQ+h.@{-x)-V^5{x) 
-(V,-iVy)d{x) 



A{Q)@{x) 


-{V, + lVy)d{x) 

-Ho-h,@{-x)-V:^b{x)^ 



(5) 



upon defining 



V2 

Ho = -—-iJ, V„ = Vo 
Zm 



(6) 



while &{x) and 5(ji:) are the Heaviside step-function and 
delta-function, respectively. Here, the barrier magnetic mo- 
ment constitutes a spin-dependent potential, where Vx = 
— pV()Cos^sin(|), Vy = — pVbsin^sin(|), — — pVocoscj). The 
intrinsic non-magnetic barrier potential is Vo, while p consti- 
tutes the effective ratio between the non-magnetic and mag- 
netic barrier, since 



P = iVl/Vo 



(7) 



where V ~ {Vx,Vy,Vz). The parameter ^ accounts for singlet 
or triplet pairing through ^ = — 1 for singlet pairing while 
^ = 1 for triplet pairing. In both cases, however, we assume 
opposite-spin pairing, corresponding to a unitary state in the 
triplet case. The exchange energy in the half-metallic ferro- 
magnet is modelled through h^, and we will later take the limit 
/zj — > fj, corresponding to a fully polarized ferromagnet. 



Solving Eq. (|5]l, we obtain the following wavefunction in 
the superconducting region: 



vsW-fJ M(es),o,o,v(es)e 



^1(7 cos Qsx 



1(7 COS 05.1: 



o,M(e,),Cv(e.)e-'^+,o 
o,Cv(7t-e5)e'^-,M(7t-es),o 



^— i^cos Q^x 



(8) 



while in the ferromagnetic region we have for an incoming 
spin-f electron with positive excitation energy e: 



^ik^ cos Qx ^ f^Q^^^^ ^-^ 



1,0,0,0 



0,1,0,0 
0,0,1,0 
0,0,0,1 



^—ik^ cos Q^x 



^ik^ cos 0^ 



gi^cose^^ 



(9) 
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In the above equations, {tg,t^} denote the transmission coef- 
ficients for electron-Uke and hole-Uke quasiparticles in the su- 
perconductor with spin a. Note that without any spin-mixing 

at the interface, one has te =1^ = in the present case. This 
is because that an incoming spin-f from the HM side can only 
be reflected normally in such a scenario without conversion at 
the interface. We comment more on this later. Moreover, q 
denotes the Fermi-level momentum in the S region while Qs 
is the propagation angle. The coherence functions are defined 
in the standard way: 



M(e) 



v(e) 







-|A(e)|2 


e 




-|A(e)P 



(10) 



We have also introduced the phase factors 

e'^± =A(e±)/|A(e±)| withe+ = e, e_ = 7:-e. (ii) 

In the normal region, we have artificially included scattering- 
coefficients for spin-l since our strategy is to do the calcula- 
tion for a strong ferromagnet, and finally take the half-metallic 
limit, corresponding to {k^,k^} — > 0. In the final expression 
for the conductance, the contribution from minority spin will 
be down a factor — > oo compared to the majority spin 
contribution, which gives us the correct result in the half- 
metallic limit. We choose q — k^, assuming that Fermi-vector 
mismatch effects simply alter the effective barrier resistance. 
In this case, 9s = 9. 

The task at hand is now to calculate the scattering coeffi- 
cients, which are needed to evaluate the conductance. To do 
so, we need to incorporate proper boundary conditions. The 
presence of a magnetic moment in the barrier, which is not 
necessarily aligned with the magnetization in the bulk HM re- 
gion, introduces new components in the boundary conditions 
as compared to the ones that mostly have been used in the 
literature. Assuming a barrier with a spin-independent poten- 
tial Vq and a spin-dependent potential of strength pVo, where 
the orientation of the magnetic moment is described by two 
angles {(|),^} as shown in Fig. [T] we may write: 

dx [^six) - \\fHM {x)] |x=o = 2mVo[l - p cos (|) (Tq (g) 03 ) 

-psin(^M(^)]VHM(0), (12) 



where we have defined the matrices 





0*(i') 



,0W 



e-'J' 
j'^J- 



(13) 



In addition, continuity of the wavefunction gives \(///m(0) — 
\|/s(0). Using these boundary conditions with Eqs. (jsj and (j9|, 
one obtains the solution for the scattering coefficients. Then, 
the conductance of the junction is expressed at zero tempera- 
ture through the dimensionless quantity 



G{eV)^ r^'d9/(9)[l + |4(ey)|2 

J-n/2 



\rUeV)\' 



(14) 



Here, /(9) is a angle-dependent factor which models the prob- 
ability distribution for incoming electrons at an angle 9. Usu- 
ally, it is chosen to /(9) = cos 9 to favor angles close to normal 
incidence, but it may also be chosen to exhibit a more discrim- 
inating tunneling cone behavior In all the plots, we will nor- 
malize the conductance on its value at voltages much larger 
than the gap, i.e. eV ^ Aq, as is usually done when compar- 
ing against experimental data since this regime corresponds to 
the normal-state conductance, and we choose /(9) = cos 9. 

Although an analytical solution for the scattering coeffi- 
cients is possible in principle, the resulting expressions are 
somewhat cumbersome, so we omit them here. For a fixed 
pairing symmetry A(9), the interface properties will deter- 
mine the behavior of the junction conductance. The inter- 
face parameters are then the spin-independent barrier strength 
Z — 2mVo I q, the ratio between the spin-independent and spin- 
dependent scattering potential p, and the orientation of the 
barrier magnetic moment determined by (|) and ^. All of these 
quantities are dimensionless. 

In what follows, we will fix the barrier strength at Z = 3 cor- 
responding to a weakly transparent interface, since this should 
correspond to a realistic experimental situation. All the inter- 
esting physics then lies in the parameters {p,(|),1'}. In the 
experimental work so far, strong sample-to-sample variations 
are seen in the results for the conductance and the critical cur- 
rent of S/HM heterostructures. As pointed out in Ref this 
is an indication that the spin-properties of the interface vary 
greatly between different samples, suggesting that the quan- 
tities {pjCj),*?} are very hard to control experimentally. The 
purpose of this paper is to obtain a fuller picture of how the 
spin-active properties of the interface influence the conduc- 
tance spectra in order to gain a clearer understanding of the 
characteristic features seen in the experimental data. 



III. RESULTS: CONDUCTANCE 

Before we proceed to a dissemination of our results, let us 
establish contact between the terminology and notation used 
in previous literature regarding proximity structures of su- 
perconductors and half-metallic ferromagnets. In our nota- 
tion, spin-mixing corresponds to the case of a finite p, but 
with (|) = 0. In this case, the effective potential felt by spin-^ 
and spin-l electrons scattered at the interface is different, and 
they pick up different phases along their scattering trajecto- 
ries. Note that the spin-dependent potential also gives rise to, 
in general, a magnetoresistance effect known as spin-filtering 
since the transmission amplitudes for opposite spins are not 
the same. Spin-flip scattering, however, requires a misalign- 
ment between the barrier magnetic moment and the bulk mag- 
netization in the half-metallic ferromagnet. In our notation, it 
is then necessary to have p 7^ and also (]) 7^ 0. Only then 
will the scattering amplitudes t\ and be non-zero in gen- 
eral, as commented on earlier From Fig. [T] it is clear that it 
suffices to vary only (|) in order to obtain both spin-mixing and 
spin-flip processes. To reduce the number of free parameters 
and still grasp the key physics, we therefore set ^ = 7t/2 in 
what follows. In this way, the barrier magnetic moment lies 
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in the y — z plane. Spin-mixing is then obtained for p 7^ and 
(|) = 0, while spin-mixing and spin-flip processes are obtained 
for p ^ and (i?^0. 

So far in the literature, the interplay between unconven- 
tional bulk superconductivity and half-metals has not been 
studied yet. We will therefore consider several bulk pairing 
symmetries in the superconducting region, including p- and 
£/-wave pairing. 



A. ^-wave pairing 

For i-wave pairing, we choose A(9) = Aq. Consider first the 
situation of pure spin-mixing, corresponding to p 7^ 0, while 
(|) = 0. This is shown in the first row of Fig. |2](left panel), 
where it is seen that the sharp coherence peak at the gap is 
replaced with broadened features upon increasing p. This is 
in agreement with Fig. 2 of Ref.^** for high values of their pa- 
rameter R, which corresponds roughly to our Z. Note that the 
subgap conductance is exactly zero, regardless of the value 
of p. The reason for this is that the usual Andreev-reflection 
where the hole has opposite spin of the incoming electron is 
not possible in the present case of a half-metal, unless spin-flip 
processes are allowed at the interface. Interestingly, Andreev- 
reflection is therefore absent in the system regardless, in fact, 
of the value of Z unless there is a magnetically inhomoge- 
neous region near the interface. We will discuss this on a mi- 



croscopic level in much more detail in Sec. Ill D 



To illustrate how spin-flip processes affect the conductance, 
we show in Fig. |2] (right panel) the case of p =0.5 for sev- 
eral values of (|). As seen, once (|) becomes non-zero, the sub- 
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FIG. 1: (Color online) The superconductor/half-metallic ferromagnet 
bilayer studied in this paper. The barrier magnetic moment may in 
general be misaligned to the bulk magnetization in the ferromagnet, 
which is assumed to be directed along z. The presence of a barrier 
magnetic moment may lead to both spin-split and spin-flip processes 
at the interface. 



gap conductance becomes finite. A large peak very close to 
eV = Ao evolves with increasing (|). However, the zero-bias 
conductance remains suppressed regardless of the orientation 
of the barrier moment. In Ref.^^, the conductance was ex- 
perimentally measured for a Pb/LaojSro.SMnOs point con- 
tact setup, where the authors found strong sample-to-sample 
variations. Some of the samples showed clear zero-bias con- 
ductance peaks, which is usually a signature of odd-frequency 
correlations or zero-energy Andreev-bound states. Other sam- 
ples displayed a clear minigap structure similar to our finding 
in Fig. |2]for the i-wave case. The strong zero-bias peak ob- 
served in the samples led the authors of Ref."*^ to speculate 
that an even-frequency p-wave bulk state was induced in half- 
metallic Lao.vSro.sMnOs by means of the proximity to super- 
conducting Pb, thus rendering the Pb/LaojSro.sMnO^ junc- 
tion into an S/S junction. Another observation in Ref."*^^ that 
supported this idea was a spectacular drop of the contact's re- 
sistance with the onset of the Pb superconductivity. 



B. p-wave pairing 

For p-wave pairing, we will consider a gap of the form 
A(9) = Aoe'^. This is known as chiral p-wave or + ipy- 
wave pairing, and is believed to be realized in Sr2Ru04. The 
gap supports Andreev-bound zero-energy states at normal in- 
cidence of incoming quasiparticles, 9 = 0, where it satisfies 
the appropriate symmetry condition A(9) = — A(7i — 9). The 
situation changes, however, in the present case where the non- 
superconducting region is half-metallic. Let us first consider 
the case of pure spin-mixing ((|) = 0) in the second row of 
Fig. [2] (left panel). As seen, the subgap conductance is still 
zero, since the gap magnitude |A(9)| is isotropic and thus pre- 
vents direct quasiparticle tunnelling into any nodes of the gap. 
Also, the usual zero-energy Andreev-bound states do not take 
part in the scattering processes since there is no possibility for 
Andreev-reflection of a spin-J, hole in the half-metal. 

The situation changes drastically once we introduce mag- 
netic inhomogeneities at the interface (right panel), corre- 
sponding to (|) 7^ 0. The subgap conductance, in particular the 
zero-bias conductance, is greatly enhanced upon increasing (|). 
The reason is that although no spin-^ holes are available in 
the half-metal, the presence of spin-flip scattering when (|) 7^ 
allows for Andreev reflection with spin-^ holes. The presence 
of Andreev reflection of majority spin holes is the reason for 
the enhancement of the conductance. 



C. rf-wave pairing 

For (/-wave pairing, we choose A(9) = Aocos(29 — 2a). For 
a = 0, there are no Andreev-bound states, while for a = 7i/4 
the order parameter supports the formation of Andreev-bound 
states in a NjtZ-wave junction. Consider first the case with 
only spin-mixing at the interface, i.e. (|) = 0. The existence 
of nodes in the gap renders the subgap conductance non-zero 
for both crystallographic orientations a = and a = 7t/4. The 
effect of increasing p is opposite for the two orientations. For 
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FIG. 2: (Color online) Plot of the normalized conductance G(eV) for a .y-wave, chiral p-wave, rf^2_,,2-wave, and d„-wa.ve symmetry in the 
rows ranging from top to bottom. In the left panels, we consider the case of pure spin-mixing for several values of p with (|) = 0. In the 
right panels, we consider additionally spin-flip processes induced by a misaligned barrier moment at the interface for several values of (|) with 
p = 0.5. 



a = 0, the conductance evolves from the typical cZ-wave bulk 
density of states profile at p = to exhibit broader features 
at p ~ 1 . For a — k/4 the conductance evolves from broad 
features at p = to a typical tZ-wave bulk density of states at 
p^l. 



Introducing spin-flip processes at the interface by allowing 
(|) ^ 0, the distinction between the two crystallographic orien- 
tations becomes clear. For a = 0, the conductance is similar 
to the bulk density of states, while for a — k/4- the zero-bias 
conductance is strongly enhanced upon increasing (|). For the 
same reason as described in the chiral p-wave case, this en- 
hancement is a result of the spin-flip induced Andreev reflec- 
tion of spin-t holes made possible by (|) 7^ 0. Interestingly, a 
large peak evolves at an energy inside the gap similarly to the 
i-wave case. One may thus ask whether the presence of a sub- 
gap peak in the conductance is indicative of surface bound- 
states induced by the spin-active interface. We shall discuss 
this question in more detail in the following section. 



D. Analytical expressions and bound-states 

In order to understand further the above results, it is instruc- 
tive to consider analytically the expression for the scattering 
coefficients. In particular, we focus on the Andreev reflection 
probability which only exists for ^-spin in the HM region. 
We find that the following general expression: 



iZ(l+pcos(^)(l+C) 

(15) 



- -4pZ/?"'sin(|) cos 9e 
X (e-'^T^ - e^'P) - cos 9(1 - C) (e"'^^ + e^'P) 
where we have defined e'P = m(9) /v(9) and 

R^z\1-pY (e"P - e-'^^)2 + 4cos4 Qe"'^ + cos^ 9 
X [e-2>AY( 1 ^ p cos + e2'P-'^T[4p (p - cos (j)) 
- ( 1 - cos^ (^)] + e'^'P (6p cos (j) + 4p2 + p^ cos^ (|) + 5) 



(16) 



Note that vanishes when either (|), p, or Z are equal to 
zero. This is physically reasonable, since the interface be- 
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FIG. 3: (Color online) Comparison of the conductance spectra for 
N|S and HM|S junctions with a li.vv-wave superconductor in the top 
row and a chiral p-wave superconductor in the bottom row. In the 
half-metallic case, we consider a spin-active interface with a misori- 
entation angle (|) = 0.5. The ratio between the magnetic and non- 
magnetic part of the barrier potential is denoted p. In all cases, the 
strong tunneling limit Z = 10 is considered. The Andreev reflection 
probability in the right panels is calculated at the peak energy in the 
dty-wave case while it is calculated at e = in the chiral p-wave case. 



comes spin-inactive in all those cases. Hence, there are no 
spin-flip processes which can mediate Andreev reflection rl. 
Consider now the triplet pairing case C = 1- I" ^ normal 
metal jchiral p-wave superconductor junction, it is well-known 
that the bound-state energies at the interface have a dispersion 
e ^ Ao sin 9. More specifically, the bound-state condition is 
given by 

2p = -Ay =71-20. (17) 

Interestingly, the Andreev-reflection coefficient Eq. ( [T5] l van- 
ishes completely for precisely these energies, regardless of the 
other parameters in the system. This is then opposite to the 
N I chiral /j-wave case where the Andreev reflection coefficient 
is unity at the bound-state energies. 

In the previous section, we pointed to the possibility that the 
emergence of strong peaks in the subgap conductance seen in 
both the s- and liw-wave cases is a signature of surface bound- 
states. This would be similar to the zero-bias conductance 
peak in N|£/„,-wave junctions originating from the existence 
of zero-energy surface states. It should be noted that an en- 
hancement of the conductance above its normal-state value 
is in general not sufficient to prove the existence of surface 
bound-states. To see this, consider e.g. a Nji-wave junction 



with a good interface contact (Z < 1), where Andreev reflec- 
tion occurs with a high probability even without any interface 
bound-states. To clarify whether an enhancement of the con- 
ductance (such as a resonant peak-structure) truly pertains to 
surface bound-states, one has to consider the tunneling limit of 
a strong barrier potential, or equivalently a low interface trans- 
parency. If the conductance is still enhanced compared to its 
normal-state value due to the presence of Andreev reflection, 
it could be a signature of resonant tunneling into a surface- 
state. We now compare the behavior of the Andreev reflection 
probability in the half-metallic limit with the corresponding 
non-magnetic case in order to acquire information about the 
origin of the conductance peak. In the top row of Fig. [3] 
we plot the conductance of a N|t/xv-wave junction (without a 
spin-active interface) and a HMjiiiy- wave junction as well as 
the respective Andreev reflection probabilities at the peak en- 
ergies. The Andreev reflection probability |r/,p is given as a 
function of the angle of incidence 9. For all plots, we have set 
Z — 10, coiTesponding to strong tunneling limit. 

In the normal metal case, the usual ZBCP is recovered and 
the Andreev reflection probability is unity for all angles of in- 
cidence. Thus, charge is transmitted into the superconductor 
as a Cooper pair via the resonant zero-energy states. Turning 
to the HM case, we wish to distinguish between the two cases 
of a weak and strong magnetic moment of the barrier. For 
p = 0.01, it is seen that the ZBCP remains, while for p — 0.5 
the ZBCP is shifted to a finite bias voltage. In the latter case, 
similar behavior was also reported for a N|ii„-wave junction 
with a spin-active interface in Ref.^^. In our case, it is nec- 
essary to have a non-zero misalignment angle (|) between the 
barrier moment and the bulk magnetization in order to gener- 
ate Andreev reflection at all, contrary to the scenario ofRef.^. 
However, the Andreev reflection coefficients shown for the 
d„-wave case in Fig. |3]indicate that the peaks cannot be as- 
cribed to resonant energy states that are available at all angles 
of incidence. In fact, the probability for Andreev reflection 
never reaches unity when p = 0.5. Still, is substantial in 
magnitude even though we are considering the tunneling limit. 
It therefore appears that surface-states are induced close to the 
interface which enable transmission processes in spite of the 
large baiTier potential, although they are not resonant in the 
sense that transmission into them occurs with a probability of 
unity. 

Turning now to the chiral /9-wave case in the lower row 
of Fig. [3] we see that the subgap conductance remains com- 
pletely suppressed when p is small. This contrasts with the 
dxy-wave case. We next increase the magnitude of the mag- 
netic part of the barrier compared to the non-magnetic part 
further, i.e. we increase p. It is now seen that the subgap 
conductance then becomes comparable to the N jchiral p-wave 
case. Note that both spin species see an effective barrier po- 
tential in the tunnelling limit even for p = 0.5 due to the large 
value of Z. Therefore, the large enhancement of the subgap 
conductance must stem from surface-induced states which de- 
cay inside the bulk. In the following section, we will employ 
a self-consistent Bogoliubov-de Gennes framework to numer- 
ically investigate whether the local DOS near the interface 
truly features such surface-bound states or not. 
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Incoming electron 



Andreev reflected hole 



, A(^) 



FIG. 4: (Color online) Andreev reflection scattering to first order for three cases: (a) N|S junction without spin-active interface, (b) HM|S 
junction without spin-active interface, (c) HM|S junction with spin-active interface. The incoming electron-like quasiparticle penetrates the 
superconducting region a distance t,s before being backscattered as a hole by the superconducting gap. 



To understand how the spin-active interface influences An- 
dreev reflection and midgap bound states on a microscopic 
level, it is useful to again compare the N|S case with the HM|S 
case. The Andreev reflection process to first order is shown in 
Fig. |4] Higher order processes may be generalized from the 
first order process along the lines of Ref.i^. From Fig. [4] 
the general recipe for Andreev reflection is seen to be a trans- 
mission of an electron-like quasiparticle into the superconduc- 
tor, which penetrates about a coherence length ^5 before it is 
scattered back as a hole by the gap A(0). Finally, the hole- 
like quasiparticle is transmitted to the non-superconducting 
region. In Fig. |4|^)' the Andreev reflection coefficient to first 

order is thus seen to be rj^'' — ti,^ijrg^i,te^e, and the contribu- 
tion from higher order processes is built along the same lines. 
We have omitted spin indices since the spin for each process is 
uniquely defined: transmission preserves spin while Andreev 
reflection flips spin. In Fig. ^h), we consider a HM | S junction 
without a spin-flip processes at the interface, corresponding 
to (|) = 0. As seen, Andreev reflection is rendered impossible 
since the spin-| hole backscattered by the gap A(6) cannot 
be transmitted into the half-metallic region due to the vanish- 
ing DOS for minority spin there. We restrict our attention to 
opposite spin-pairing superconductors only, such as .s-wave, 
chiral p-wa\e (S^ — 0), and li-wave. For an equal spin-pairing 
superconductor, Andreev reflection is obviously possible even 
without any spin-flip processes at the interface. In Fig. |4|c), 
we consider a HM|S junction with spin-flip processes at the 



interface. As mentioned previously, Andreev reflection is now 
possible even for a backscattered minority spin hole due to the 
spin-flip probability at the interface. 

The microscopic picture shown in Fig. |4] also allows us 
to understand how the midgap bound states are influenced 
by the spin-active interface. To do so, we first briefly reca- 
pitulate the results of Ref for the N|S case without spin- 
active processes. Focusing on £ = where the retroreflec- 
tion property holds in the S region, one can calculate the to- 
tal probability for Andreev reflection obtained by summing 
all orders of scattering diagrams such as the ones shown in 
Fig. |4] In doing so, the total probability is proportional 



to I ir=ok 



2n\ 



]" I , where j± represents the phase 



contribution from the internal phase of the superconducting 
order parameter while \r\ = \re^e\ — k/i^/il- For a dxy-wave 
superconductor where e''^-^^+' = —1, it is seen that all orders 
sum in a coherent way and the total Andreev reflection prob- 
ability can be shown to equal unity. In the present case, the 
phases picked up by the scattered particles are longer spin- 
degenerate. In particular, we see from Fig. |4|c) that the prob- 
ability for Andreev reflection to first order is equal to 



„(i) 



(18) 



The crucial point is now that, whereas the branch-converting 
reflection coefficients r^^h have spin-independent scalar 
phases without a spin-active interface, they are spin- 
dependent otherwise. In the former case, one has rj.o^A,-o ~ 
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g-iy(e)+ii> jjjj^j ri,„^^ _„ ^ e'T(^)-'* with a =t,i. When sum- 
ming the Andreev reflection processes to all orders, one ob- 
tains products of rea^h-o and r/,(j^f.-a which effectively 
gives a phase-factor e'(^-^''+^ while the other scalar phases 
cancel each other. When the interface is spin-active, the scalar 
phases become spin-dependent, and one effectively gets an 
additional contribution AiS^o = "^o — iS^-o in the phase of the 
effective Andreev reflection coefficient. For this reason, the 
summation over all orders n is altered and the resonant states 
at e.g. e = in the dxy-v^s.Mt case are shifted. The spin- 
dependent phase-shifts ■do depend on both p and (|) and are 
the reason for why the conductance is qualitatively altered in 
the presence of a spin-active interface as shown in both Fig. [2] 
and Fig. |3] 



IV. PROXIMITY EFFECT AND LOCAL DENSITY OF 
STATES: A LATTICE STUDY 

In this section, within a self-consistent scheme of compu- 
tation, we investigate the proximity effect and its influence 
on both the superconducting order parameter and the local 
density of states in the HM|S bilayer for a model system as 
described in the Fig. |5] The analysis is based on the case 
of a junction configuration with the barrier aligned along the 
y direction in the xy plane. The pairing amplitudes are ex- 
pressed in terms of the components along the x and y axes of 
the square lattice. In particular, within such configuration, due 
to the symmetry properties of the examined order parameters, 
we do expect that Andreev bound states are effective only for 
the chiral p-wave paired state (for the analysis on the lattice 
we do not consider the case of d^y symmetry of the supercon- 
ducting order parameter). The formalism to be used when cal- 
culating the order parameter profile and the density of states 
is based on a lattice BdG-approach similar to that adopted in 
but extended to the case of a spin active interface. The 
total Hamiltonian H of the system may be written as 



H — Hp + Hs + Hj + Hi, 



(19) 




FIG. 5: (Color online) Schematic description of the lattice geome- 
try for the HM|S bilayer junction indicating the reference axis sys- 
tem, the notation for the size (2Lt x 2Lv), the position of the barrier 
{x = 0), the sketch of the pairing configurations (d-wave, p-wave and 
local s-wave), respectively. On the left side, we show a sketch of the 
density of states for the half-metal. 



where Hf and Hs accounts for the ferromagnetic and super- 
conducting layers, while Hp and Hi describe the tunnelling 
processes and the scattering potential at the interface region. 
We have 



Hi 



- L(i,j> ,o ha (CjoCjo + h.c.) + ?7A«itnii 

- I(i,j> Vk (niinji + "iinjt) - A* Li,o «ia 



(20) 



where (i,j) denotes nearest-neighbor sites, {c-^jQa} are cre- 
ation and annihilation operators of an electron with spin a on 
site i = {ix,iy), while is the number operator. We take the 
exchange field to be non-zero only on the ferromagnetic side, 
where we let hf — > /j in correspondence with our assumption 
of a half-metallic limit. The hopping amplitudes are chosen 
such that I = fF|=fs| — tsi=t. Above, /j is the chemical po- 
tential, while Ui and V/ denote on-site and nearest-neighbor 
interaction on side I. 

The two layers communicate by means of the tunneling 
Hamiltonian term, which reads 



Hj 



h.c. 



(21) 



where the sites (i, j) are located at the surface of the F and S 
layer Finally, the scattering potential at the interface is mod- 
elled through the term 



Hi 



iap 



(22) 



where {a, p} are spin-indices, while i denotes a lattice site 
on the surface of the F or S layer (for convenience the inter- 
face has been placed at the site x = as indicated in Fig. [5]l. 
Here, Vo is a spin-independent scattering potential, roughly 
corresponding to the parameter Z introduced previously, while 
Ym is a spin-dependent scattering potential which gives rise to 
spin-mixing and spin-flip processes. The Pauli-vector matrix 
is given as o = {ax,Oy,0^) and 1 is the 2x2 identity ma- 
trix. To facilitate comparison with the notation and parameters 
used when calculating the conductance, we similarly define 



Vm = — pVo(cos^sin(|), sin1'sin(|),cos(|)), 



(23) 



such that p denotes the relative weight of the spin-independent 
and spin-dependent potential while {(j), ^} provides the direc- 
tion of the magnetic moment at the interface. We analyze 
the effects of the spin-active barrier on the amplitude and the 
phase of the superconducting order parameter by solving the 
BdG equations on the lattice within the extended Hubbard 
model introduced above. The calculation is performed for 
the case of a planar bilayer junction of dimension IL^xlLy. 
We have considered lattice sizes of Lx=Ly= 40, 50, 60 (lattice 
constant is the the unit of length). The results obtained do 
not show qualitative nor significant quantitative changes for 
these values of Lx,Ly. The case stiudied here corresponds to 
superconducting and magnetic coherence length of the order 
of ^5 ~ 7 and ^//^ — 1 • Hence, size effects can be considered 
negligible for the systems we consider Concerning the ratio 
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i^s/'^HM, Other computations have been performed at differ- 
ent values of the pairing strength V , thus varying ^5. These 
cases do not show quahtative changes in the results. In this 
framework, the modification of the coherence length is lim- 
ited by two conditions: i) the requirement of stable supercon- 
ducting solutions in the phase diagram for the order param- 
eters in the desired symmetry depends on the pair coupling, 
and ii) the computational demand is related to the size of the 
matrix Hamiltonian. Hereafter, the discussion will focus on 
the case L^=40. From Eq. |20]and the BdG formalism already 
described in Ref.^^, but extended to the case of a spin active 
barrier, we compute the spatial variation of the superconduct- 
ing order parameter for different pairing symmetry and as a 
function of the barrier parameters. Due to the presence of a 
spin active barrier one has to introduce a four component Bo- 
goliubov basis on each atomic site to take into account both 
particle-hole spin flip processes as well as the pairing channel 
of particle-hole resonance. This introduces an extra factor in 
the computational complexity. In particular, among the var- 
ious results obtained, the focus is on the modification, due 
to the split exchange and the spin flip coupling, of the su- 
perconducting order parameter evaluated at the HM — S inter- 
face versus ((|),Vm,Vo)- In doing that, we have to properly 
choose the interaction strength both in the magnetic as well 
as in the superconducting subsystem of the junction in order 
to get the desired microscopic quantum states. To this end, 
the effective exchange amplitude in the ferromagnetic region 
is taken as hp = 4. Of to have a half-metallic behavior and a 
profile for the z-component of the magnetization with zero 
spin minority carriers. Furthermore, to get an i-wave, li-wave 
and a chiral p-wave symmetry within the superconducting re- 
gioon, three different sets of attractive pairing amplitudes and 
chemical potentials have to be considered^. For the onsite 
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FIG. 6: (Color online) Color map of the s-wave on-site pairing am- 
plitude at the interface As[0] with respect to the bulk value g as 
a function of the misalignment angle (|) and the magnetic barrier 
strength Vm- The non magnetic scattering potential has an amplitude 
Vo = 2t. 



s-wave, we assume a value of the chemical potential equal to 
/J = -0.2f , with Uf=0,Vf=0 for the F side and t/s = - 1 .5f , 
Vg — for the S side. For the chiral /7-wave, we choose the 
value /J = — 1.6f, with Up ^ 0, Vp = within the F side and 
Us = 0,Vs = —2.5t for the S side. For the cZ-wave, we choose 
Ij = -0.2f , with Up =0,Vf =0 within the F side and Us = 0, 
Vs = -2.5f in the S side. 

Furthermore, the tunnelling matrix element is kept fixed 
and chosen equal to tp = t. 

The interaction terms in Hp and Hs are decoupled by means 
of a standard Hartree-Fock approximation such that the mag- 
netic and pairing channels originate from the on-site and the 
intersite interactions, respectively: 



Up [(ni;)nit + («it>nii - ("it) («i;>] 



Us 



Ai44 + A?Ci|Ci|-|Ai| 



Here, we have introduced the on-site Aj = (ci|Ci|) and the 
bond pairing amplitude on a bond Ay = (ci|Cj|), with the av- 
erage {K) yielding the expectation value of the operator K 
over the ground state. Moreover, the on site z-component 
OT,i = 5((ni|) — {nil)) (ji:,37)-components of the mag- 

netization m^i — j((c[|Ci| +/i.c.)),mvi = j((c[|Ci| — /i.e.)) are 
iteratively determined up to the required accuracy to get the 
spatial dependence of the vector spin polarization. From the 
pairing amplitudes, it is possible to construct the supercon- 
ducting profiles for the different symmetries {d- and p-wave) 
in the singlet (S) and triplet (T) channel in terms of compo- 
nents of the z-projected axial spin operator. They are defined 



A./(i) 
Ap.(i) 
Ap.(i) 



= A 



i,i+x 



-A|t,-A<£,-A<t,)/4 



(a!2,-a|i!,)/2 
(a!2,-a2,)/2 



for d^2_^2, Px and p^-wave, respectively. Here one has to in- 
troduce the singlet and triplet pairing amplitudes on a bond, 
given by 

Afj = (Aij+Aji)/2 
Aj - (Aij-Aji)/2 . 

We adopt open (periodic) boundary conditions for the direc- 
tion X perpendicular {y parallel) to the interface, taking the 
Fourier transform due to the translational invariance along the 
y direction of the relevant physical quantities. 

In order to understand the role of the spin active barrier in 
tuning amplitude and phase of the different symmetry order 
parameters we have investigated their evolution in terms of 
the barrier parameters at zero temperature. The angle (|) tunes 
the direction of the barrier spin moment from parallel to the 
z-quantization axis at (|) = to being parallel to the x- direction 
for (|) = 7t/2. Different values of ^ do not change the results 
due to the x — y symmetry in the spin space. Thus, it suffices 
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to consider the case ^ = 0. Furthermore, we have chosen 
one representative case for the regime of nonmagnetic barrier 
strength that corresponds to a situation of reduced electron 
density at the interface compared to the electron distribution 
in the HM and S sides. In this respect, the way of introducing 
a scattering potential at the interface may lead to extra effects 
if compared to the delta potential considered in the BTK for- 
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FIG. 7: (Color online) Color map of the real part of the p v component 
of the pairing amplitude evaluated at the interface A^^ [0] with respect 
to the bulk value Ap^ s as a function of the scaled misalignment angle 
and the magnetic barrier strength Vm/Vq. The non magnetic 
scattering potential has a given amplitude of Vq = 2t. 
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FIG. 8: ((Color online) Color map of the imaginary part of the py 
component for the pairing amplitude evaluated at the interface Ap^, [0] 
with respect to the bulk value Ap^ .g as a function of the scaled mis- 
alignment angle 2<^/n and the ratio of the magnetic barrier strength 
with respect to the nonmagnetic one, Vm/Vq. The non magnetic scat- 
tering potential has an amplitude Vq = 2t. 



malism, since the amplitude of the non magnetic potential de- 
termines the average electron occupation at the barrier site as 
well as in its proximity. 

s-wave pairing As one can see in Fig. [6j the value of the 
s-wave order parameter does not vary significantly in the full 
range of values for Vm and (|). Thus, the ampUtude reduction 
with respect to the bulk value is basically controlled by the 
presence of the half-metallic ferromagnet. Only by approach- 
ing the regime of Vm — Vo the effects of the spin active barrier 
become more relevant leading to a strong suppression of the 
pairing amplitude. In this case it is possible to distinguish two 
different behaviors corresponding to the spin-mixing or spin- 
flip barrier regime. Spin-mixing effects (i.e. (|) = 0) are not 
much relevant for the s-wave proximity effect as the pair am- 
plitude exhibits only a slight reduction as one tunes the spin- 
dependent scattering from the regime p = to p ~ 1 . This can 
be understood because the change of Vm tends to reinforce the 
magnetization even at the barrier site. There, the proximity 
between the half-metal ferromagnet and the superconductor 
leads to a matching of the magnetization from full- to zero- 
spin polarization in moving through the interface. Otherwise, 
spin-flip mechanisms lead to a larger reduction of the pair am- 
plitude when approaching the limit p — > 1. In this case, the 
increase of the transverse magnetization (parallel to x in the 
spin space) at the barrier site leads to extra scattering for the 
singlet pairs that in turn sums up to the pair breaking effect 
due to the presence of the half-metal ferromagnet in a way to 
get about a 70% reduction. 

Chiral p-wave pairing Consider next the case of chiral p- 
wave pairing (i.e. A ^ p^ + ip^) on the superconducting side 
of the junction. The chiral state exhibits time reversal sym- 
metry breaking for a spin triplet configuration. It is well 
known that the interface properties can be quite unusual even 
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FIG. 9: (Color online) Color map the d^2_^a pairing amplitude at the 
interface A^ [0] with respect to the bulk value A^ g as a function of 
the scaled misalignment angle 2^/ti and the ratio of the magnetic 
barrier strength with respect to the non-magnetic one, Vm/Vq. The 
non magnetic scattering potential has an amplitude Vo = 2t. 
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at the boundary with the vacuum due the possibility of emer- 
gent exotic edge states. Here, we analyze the consequences 
of spin-mixing and spin-flip barrier processes on the two p- 
wave components separately. Though the gap amplitude is 
isotropic in fc-space, as for the i-wave symmetry, the pairing 
amplitudes along x and y evaluated at the interface exhibit a 
completely different behavior. There are many distinguishing 
features that can be extracted from inspecting Figs. |7]and[8] 
Regarding effects induced by spin-mixing and spin-flip pro- 
cesses, note that i) spin-mixing and spin-flip reduce the pair 
breaking effects due to the proximity with the half-metal, ii) 
the reduction is not equivalent for the x and y components, 
iii) the pair breaking effects are more pronounced for the px 
than the py component. Indeed, for the case i) one notices 
in Figs, [t] and [S] that in the regime Vm/Vq ^ 1, the pair am- 
plitude tends to grow at any given angle with a slope that is 
more pronounced for the p^ component as compared to the py 
one. Concerning the point ii), the maximum of the pairing am- 
plitude occurs at the phase diagram positions individuated by 
(c^, Vm/Vo) - (0, 1) and {^,Vm/Vq) = (7i/2, 1) for the p^ and 
the Py amplitude, respectively. We argue that the presence of 
Andreev bound states in the spectrum, due to the change of 
sign of the px component in the direction perpendicular to the 
interface, leads to a more significant barrier influence of the 
correspondent pairing component with respect to the py one. 

d-wave pairing In Fig. |9] we show the evolution of the 
interface d-wave pairing amplitude in terms of the scaled an- 
gle and scattering barrier parameters. Recall that for the cho- 
sen junction geometry there are no Andreev bound states in 
proximity to the half-metal. Hence, even though the pairing 
is anisotropic, the absence of Andreev states leads us to ex- 
pect a behavior similar to the i-wave case. Indeed, the profile 
of the c/-wave pairing amplitude at the interface is similar to 
that of the isotropic i-wave with an increase of the pair break- 
ing effects in the regime of large spin-mixing and spin-flip. 
However, it is worth pointing out that close to the regime of 
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FIG. 10: (Color online) Density of states for the chiral p-wave eval- 
uated at one representative site position (; = 2) within the range of 
one superconducting coherence length from the barrier (i = 0). The 
non magnetic potential has a value of Vq = 0.5t. 



maximal spin-flip scattering ((|) = 7t/2), where the orientation 
of the magnetic moment at the interface is perpendicular to 
the easy axis of pair formation, the combination of the spin- 
flip and the presence of the half-metal leads to a global phase 
change in the d-wave amplitude. This is reminiscent of the 
oscillating behaviour one would expect in a conventional fer- 
romagnet/superconductor junction. Here, it is the anisotropy 
of the pairing that favors the formation of small oscillations in 
proximity of the barrier. 

Interface density of states (DOS) In the previous section, 
we saw how a finite subgap conductance was obtained even 
in the strong tunnelling limit when the superconducting pair- 
ing was unconventional and the interface allowed for spin-flip 
processes. We speculated that the physical explanation behind 
this phenomenon was the generation of surface bound-states, 
appearing due to an interplay between the spin-active proper- 
ties of the interface and the internal phase of the superconduct- 
ing order parameter In the absence of a spin-active interface, 
we found that subgap tunnelling vanished completely. We 
now investigate whether such surface bound-states are truly 
present or not within the lattice BdG-model, by focusing on 
the p-wave chiral type of pairing. This is expected to yield 
Andreev bound states for the geometry given in Fig|5] 

In Fig 10 we show the behavior of the total DOS 



evaluated at one representative site position within a range of 
one superconducting coherence length from the barrier. There 
is a substantial difference in the results for the total DOS ob- 
tained when comparing the case of a complete non-magnetic 
barrier ((|) = and Vm — 0) with the case of a pure spin-flip 
barrier scattering potential ((|) — n/2 and a varying Vm)- The 
energy has been rescaled with respect to the gap amplitude 
/ig determined inside the superconductor within the same 
formalism, that is at a position where the DOS exhibits a full 
gap as expected for the chiral p-wave symmetry and the order 
parameter is uniform and unaffected by the interfaces. Note 
that for non-zero (|) the change in the amplitude of Vm leads to 
extra midgap edge states. The energy of these states depend 
strongly on the magnetic barrier potential. Indeed, they ap- 
pear at the edge of the gap for small values of Vm and as the 
amplitude of the magnetic moment at the barrier increases by 
tuning Vm, they shift towards low energies. The asymmetry of 
the resonant spectra is related to the presence of a small, but 
finite spin polarization proximate to the barrier on the super- 
conducting side. 



V. SUMMARY 

In summary, by means of continuum and lattice BdG for- 
malisms we have computed the conductance spectra, super- 
conducting order parameter, and the density of states of a bi- 
layer system made of a half-metal and an unconventional su- 
perconductor which are brought into contact through a spin 
active tunneling. These quantities have been computed for dif- 
ferent forms of the pairing amplitude both for spin singlet s- 
and (i-wave as well as spin triplet chiral p-wave symmetry. We 
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have shown how the spin-flip and spin-mixing scattering pro- 
cesses at the interface influence the structure of the supercon- 
ducting order parameter in the case of different pairing sym- 
metries, namely s-wave, chiral p-wave, and J-wave pairing. 
These scattering prcoesses lead to different charge transport 
features such as subgap conductance and midgap Andreev as- 
sisted tunnelling due to resonant states, depending on which 
pairing symmetries are considered. The novel subgap fea- 
tures we find are present only for a non-zero misalignment be- 
tween the half-metalhc magnetization and the magnetic mo- 
ment of the barrier. The energy position of the Andreev as- 
sisted charge processes turns out to be sensitive to changes in 
the ratio between the magnetic and non-magnetic scattering 
potential at the interface. We have also studied the pairing am- 
plitude in a lattice model and computed its dependence at the 
interface on the barrier scattering strength, as well as its de- 
pendence on the misalignment angle between the half-metal 
magnetization and the barrier spin moment. For the case of a 
spin triplet chiral /j-wave, midgap states at non zero misalign- 
ment angle between the half-metal magnetization the barrier 
spin moment have been found by computing the density of 



states of the superconductor close to the interface. A hallmark 
of such edge states is the strong dependence of the ratio be- 
tween the magnetic and non-magnetic scattering potentiiil in 
the subgap DOS. Finally, we point out that the present study 
reveals highly non-trivial features of the proximity effect be- 
tween a half -metal and an unconventional superconductor in 
the presence of a spin active interface even without invoking 
the occurrence of exotic mixed parity pair components. The 
present analysis is the starting point for further investigation 
of the role played by induced or subdominant pairing ampli- 
tudes in heterostructures based on half-metal and imconven- 
tional superconductors. 
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